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We compute the three-point correlation function of primordial scalar density perturbations in a 
general single-field inflationary scenario, where a scalar fleld has a direct coupling with the Ricci 
scalar R and the Gauss-Bonnet term Q. Our analysis also covers the models in which the Lagrangian 
includes a function non-linear in the field kinetic energy X — —{d4>)^/2, and a Galileon-type field 
' ^ , self-interaction G{(j>, X)0(j}, where G is a function of tp and X. We provide a general analytic formula 

I ' for the equilateral non-Gaussianity parameter /^l"'' associated with the bispectrum of curvature 

I perturbations. A quasi de Sitter approximation in terms of slow-variation parameters allows us to 

. derive a simplified form of /^l"'' convenient to constrain various infiation models observationally. 

If the propagation speed of the scalar perturbations is much smaller than the speed of light, the 
Gauss-Bonnet term as well as the Galileon-type field self-interaction can give rise to large non- 
. Gaussianities testable in future observations. We also show that, in Brans-Dicke theory with a 

' field potential (including f{R) gravity), /^l"'' is of the order of slow-roll parameters as in standard 

\^ ' infiation driven by a minimally coupled scalar field. 

I ! 

PACS numbers: 98.80.Cq, 95.30.Cq, 04.62. +v 

O: 

^ ■ I. INTRODUCTION 

(-^ . 

O ' The idea of cosmic acceleration in the early Universe P, 01 was originally introduced as a way of addressing 
I , the flatness and horizon problems plagued in the big bang cosmology. The inflationary paradigm can provide a 
2 ■ causal mechanism for generating density perturbations responsible for the Cosmic Microwave Background (CMB) 
' temperature anisotropics Q. The standard, single- field slow- roll models of inflation predict adiabatic and Gaussian 
i primordial perturbations with a nearly scale-invariant spectrum. This property is consistent with the observed CAIB 
I— I ' anisotropics 0-01 as well as the large-scale structure data |8i] . 

Over the past 30 years, many models of inflation have been proposed in the framework of particle physics or extended 
^ ■ theories of gravity (see e.g., [9| for reviews). From the information of the spectral index n-ji of the scalar perturbations 
', as well as the tensor-to-scalar ratio r, the CMB observations by the WMAP satellite have been able to rule out some 
of those models [E-0| • While the observables n-jz and r are derived in linear cosmological perturbation theory, it is also 
possible to distinguish between a host of inflationary models further by comparing the non-Gaussianity of primordial 
perturbations with the WMAP data [lol - [l6j . In particular it is expected that the PLANCK satellite 17] will bring 
■ us more precise data of non-Gaussianities within a few years. 

, The amount of non-Gaussianity can be quantified by evaluating the bispectrum of curvature perturbations TZ, as 

(7^(fcl)7^(fe2)7^(fc3)) - (2^)35(3)(fel+fc2 + fc3)B7^(fcl,fc2,fc3), (1) 

where TZ{k) is a Fourier component of TZ with a wave number k. Conventionally the bispectrum B-ji translates into 
' a non-linear parameter /nl in order to confront theoretically predicted non-Gaussianities with observations [13] ■ 
H There are two different shapes of the primordial bispectrum: (1) "local" type [H, El], and (2) "equilateral" 
' type jl^l- The first one arises from a local, point-like non-Gaussianity given by TZ{x) = TZh{x) -\- (3/5)/^f 

where TZi^ is a linear Gaussian perturbation. In single-field inflation models driven by a sl owly varying potential, the 
predicted value of /^l^' is of the order of slow- roll parameters, i.e. I/nl'^'I ^ 1 [Hl-O IH, HJ- This small amount of 
non-Gaussianity is consistent with the WMAP 7-year bound: -10 < f^l'^^ < 74 (95 % CL)^. 

In the equilateral shape of non-Gaussianities, the momentum dependence of the function B-jijki, fc2, fca) arising 
from non-canonical kinetic term models ( "k- inflation" |22l|) can be approximated in a suitable form |20j . In k-inflation 
models, including Dirac-Born- Infold (DBI) inflation f23'| and (dilatonic) ghost condensate [23|, one can realize the non- 
linear parameter /^l"'' larger than the order of unity provided that the propagation speed Cs of scalar perturbations is 
much smaller than 1 (in unit of speed of light) The WMAP 7-year bound is -214 < /"f < 266 (95 % CL) 

@, but it is expected that the PLANCK satellite can reduce this limit by one order of magnitude. This will provide 
us an opportunity to distinguish k-inflation from standard slow-roll inflation. We also note that it is possible to give 
rise to large non-Gaussianities in multiple scalar-fleld models [31 
models [ia, |47| , models having a preheating stage after inflation 
m (see also Refs. [50.])- 



[401 ■ curvaton models [4lH45j . modulated reheating 
48| , and models with a temporal non-slow roll stage 
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In this paper we shall compute the bispectrum B-jz and the equilateral-type non-linear parameter /^l"' for general 
single-field inflation models described by the action ([5]) below. In low-energy effective string theory there is a scalar 
field (j) called dilaton coupled to the Ricci scalar R with a form F{(j))R [5l|. A field coupling of the form £,{(t>)g, where 
Q is the Gauss-Bonnet term, also arises as a higher-order string correction to the low-energy effective string action 
[5^ . Furthermore the higher-order string correction contains a non-canonical kinetic term like (dtj))^ as well as a field 
self-interaction of the form g{(p){d(p)'^Dip in the action. For constant g{(j)) the latter is linked to the Lagrangian of a 
covariant Galileon field that respects the Galilean symmetry 9^0 d^cj) -\-b^ in the Minkowski space-time [s^ 
The cosmological dynamics in the presence of the Galileon-type interaction g{(/)){d(f>)'^0(f> have been extensively studied 
recently in the context of inflation [5^ [s^ and dark energy [s^l ■ We accommodate non-linear field derivative terms 
as the Lagrangian of the form P((/), X) — G{(j), X)n0, where P and G are functions in terms of </> and X = —(90)^/2. 
Note that non-Gaussianities in the models described by the Lagrangian F{<j))R + P{(j),X) without the terms ^{(l>)G 
and G{(j),X)D<j) were recently studied in Ref. [ssj . 

Our action ([2]) can be also viewed as describing a kind of modifled gravitational theories [11]. In fact, this covers 
the so-called scalar-tensor theories [13] such as Brans-Dicke (BD) theory (Hlj. Since metric f{R) gravity is equivalent 
to BD theory with the BD parameter wbd = ^H, our results can be applied also to the Starobinsky's inflation 
model f{R) = R + aR^ Q. Since the field propagation speed Cs is unity in BD theory, we will show that I/nl"' I ^ 1 
as in conventional slow-roll inflation. On the other hand, the terms £,{(j>)G and X)D(t) as well as P(0, X) can give 
rise to Cs much smaller than 1. In such theories it is possible to realize large non-Gaussianities detectable in future 
high-precision observations. 

This paper is organized as follows. In Sec. [Tl] the field equations of motion are derived for the action ([2]) on the flat 
Friedmann-Lemaitre-Robertson- Walker (FLRW) background. In Sec. IIIII we obtain the second-order action for the 
curvature perturbation TZ and present the solution of its mode function at linear level on the de Sitter background. 
In Sec. lIVI the third-order perturbed action is derived in order to compute the three-point correlation function of 7?. in 
the interacting Hamiltonian picture. We also present a general analytic formula of /^l"'' valid in the quasi de Sitter 
background. In Sec.|V]we provide a simpler expression of /^l"'' under the expansion of "slow- variation" parameters. 
This approximate formula is enough to estimate the amount of non-Gaussianities in practical purpose. In Sec. IVII we 
apply our general results to a number of concrete models of inflation: (1) k-inflation, (2) generalized Galileon model, 
and (3) Brans-Dicke theory. Sec. I VIII is devoted to conclusions. In Appendix we present the detailed procedure to 
derive and manipulate the third-order perturbed action. 



II. THE MODEL AND BACKGROUND EQUATIONS 

We start with the following action 



S = d x^—g 



-^F{^)R + p(0, X) ~ mg ~ G(0, x)n0 



(2) 



where g is the determinant of the space-time metric 5^1/, = (SttG)^^/'^ is the reduced Planck mass (G is gravi- 
tational constant), and is a scalar field with a kinetic term X = —{l/2)g^'^d^(j)dv(j)- The functions F{(f)) and ^((/>) 
depend on (f) only, whereas P{(f>,X) and G{(j),X) are functions of both </> and X. The field couples to the Ricci 
scalar R as well as the Gauss-Bonnet term Q defined by 

g = R^- 4R^p R"^ + Ro^p^sR"^''^ , (3) 

where Rap is the Ricci scalar and Rap-yS is the Riemann tensor for the metric g^^^. As we mentioned in the Introduction, 
the action ([2]) covers a wide variety of single-field inflationary models. 

We consider the flat FLRW space-time with a scale factor a{t), where t is the cosmic time. Then the background 
equations are given by 

El = ZMlFH^ + UllHF + P- 2XP^x ~ 2AH^i - 6H^XG,x + 2AG,0 = , (4) 



E2 = Mll^FH^ + 2Ml^HF + 2MlyFH + MI^F + P- WH^i - IQHHi - 8H^i - G,x<i>^4> - G,^^^ = , (5) 



= (Fx + 2XPxx + 6ff<^G,x + GH^XG^xx - 2AG,^x - 2G,^)0 

^(f,x — 6HG,^ — G^, 

where H = d/a is the Hubble parameter, and a dot represents a derivative with respect to t. These equations are not 
independent because of the Bianchi identities, i.e. (pE^ + Ei+ 3H{Ei — E2) = 0- In Eq. ([5]) the term G,x4>'^4' vanishes 



P^ - 6M^,H^F^^ - iMlfiF^^ + 2mX^ + 2m^H£,^^ = , (6) 
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for the theories with G x = 0, in which case the Lagrangian G(0)n0 can be regarded as a part of the Lagrangian 
P{(f>,X). If G X 7^ 0, however, the term G{(l), X)D(j) should be treated separately from the term P{(j>,X). 

In order to derive the slowly varying parameter —H/H^ we consider the combination {E2 — Ei)/{Mp^H^F) = 0, 
which gives 



H F F XPx 4:H(, 8H^ 4^ 3(I)XG x (t>XG,x 2XG , 



2HF 2H^F M^^H^F M^^F M^^HF M^^F M^^HF M^^H^F M^^H^F 



(7) 



During inflation the Hubble parameter changes slowly, so that the condition e <C 1 is satisfied. Hence, in general, 
we require that each term on the r.h.s. of Eq. (O is much smaller than unity. In Sec. |V] we shall use this property 
to obtain a simple expression for the equilateral non-linear parameter /^l"' ■ ^'^^ later convenience we introduce the 
following "slow-variation" parameters 

III. SECOND-ORDER ACTION AND LINEAR PERTURBATIONS 

In order to compute primordial scalar non-Gaussianities we need to expand the action 1^ up to third order in the 
perturbations, taking into account, up to a gauge choice, both the perturbations in the scalar field, 5(j), and in the 
scalar modes of the metric. The interacting Hamiltonian of perturbations follows from the third-order Lagrangian, 
by which the three-point correlation function of curvature perturbations TZ can be evaluated in the framework of 
quantum field theory (T^j- In order to calculate the vacuum expectation value of the correlation function, the mode 
function of TZ should be known in the quasi de Sitter background. Once the mode function is obtained at linear level, 
one can derive the power spectrum of TZ generated during inflation. The linear perturbation equation for TZ is known 
by the second-order perturbed action. 

For the derivation of the action expanded up to third order in the perturbations, it is convenient to work in the 
ADM formalism ^] with the line element 

ds^ = -N^dt^ + h,j{N'dt + dx'){N^dt + dx') , (9) 

where N and N"^ are the lapse and shift functions, respectively. Here we consider only scalar metric perturbations 
about the flat FLRW background. In doing so we expand the lapse N and the shift vector as iV = 1 -I- a and 
Ni = diip, respectively {i = 1,2,3). In fact, in the ADM formalism, the lapse N and the shift N'^ are Lagrange 
multipliers, so that it is sufficient to know N and A^' up to first order. This is because the third-order and the 
second-order terms in N, iV* multiply the constraint equations at zero-th order and at first order, respectively, so that 
their contributions vanish (12.. 66] . We choose the uniform-field gauge with S(j) = 0, which fixes the time-component 
of a gauge-transformation vector We gauge away a field E that appears as a form E^ij inside /ly , by fixing the 
spatial part of . Then the three-dimensional metric can be written {t)e'^''^Sij. This results in the following 

metric 

ds^ = - [(1 + a)2 - a-2(f)e-2^(aV')2] dt^ + 2^^t^:dtdx' + a^{t)e^'^{dx^ + dy^ + dz^) , (10) 
se SI 

theory, that is [6J] 

ds^ = -{\ + 2a)df ^2d^ipdtdx^ + a^{t) {1 + 2TZ) {dx^ + dy^ + dz^) . (11) 
Expanding the action ^ up to second order, one finds 



where (dtp)'^ — (dtpi) {d4>i) = {dx'4')'^ + (dytp)'^ + {dztp)'^ ■ Here and in the following, same lower latin indices are summed, 
unless otherwise specified. At linear level the metric (|10p reduces to the standard one used in linear perturbation 



dt d'^x a-^ 



■ 2w^ ■ „9 Wn „T 2wi „o • 1 9 U'4 „ „ 

-iwiTZ^ + -^TZd^ij; ^aS^ ^ad^TZ + iw2aTZ+ -ws.a^ + diTZ diTZ 

3 



(12) 



where 



= M^iF~8H(, (13) 

W2 = Ml^{2HF + F)-2^XG,x-2AH'^i, (14) 
W3 = ^QM^^FH^ - QM^^HF + 3(APx + 2X^P^xx) + UAH^i 

+ l8H<j){2XG^x + X^G^xx) - 6{XG.^ + X^G ,<t>x) , (15) 

= M^,F~8^. (16) 



4 



In the action ([12]), both the coefficients of the terms aTZ and TZ^ vanish by using the background equations of motion. 
This is related to the fact that the field TZ does not have an explicit mass term. Furthermore, in ([T^. the term 
quadratic in ip vanishes after integrations by parts. The equations of motion for and a, derived from (|12[) . lead to 
the following two constraints 



2W3 
3^2 



a + 37^- 



2wi 1 „2 



W2 



where 



2(Ml,F - 8Hi) 



2wi _ 



(17) 
(18) 



(19) 



In k-inflation with F = 1, ^ = 0, and G = one has Li ^ 1/H. In general cases, expansion in terms of the 
slow- variation parameters defined in Eq. ([S]) gives 



Li 



1 

H 



1-Isf + 4(5^ + Sgx + O(e') 



Plugging the relation a = (2wi/u>2)7t into Eq. (IT2|) and integrating the term TZd^TZ by parts as 

c{t)lZd'^n = c {dnf/2 + total derivatives , 

one finds 



5*2 = / dtd-'xa^Q 



where 



Wi{4wiW3 + 9W2) 



3wj 



3{2wlw2H — W2W4 + AwiWiW2 ~ 2W1W2) 
'Wi{4:WiW3 + 9^2) 



(20) 

(21) 
(22) 

(23) 
(24) 



More explicit expressions for Q and are given in Appendix [Xj In order to avoid the appearance of ghosts and 
Laplacian instabilities we require that Q > and > 0, respectively. Using Eqs. pT)) and (^5]) . we can rewrite 
Eq. (dll) as 



ip = -LiU + x, where a^x = a^— ^• 
We introduce a parameter defined by 

_ Qc^ 2W1W2H — W2W4 + 4:WiWiW2 — 2w1w2 



(25) 



(26) 



In k-infiation with F = 1, ^ = 0, and G — i) this reduces to the slow-roll parameter e = —H/H^. In general, one can 
expand ts in terms of slow- variation parameters, as 



cs = e + ^<5f + 5gx - 4^5 + 0{e^) . 



(27) 



The equation of motion for TZ follows by varying the Lagrangian £2 — a^Q[7t^ — (c^/a^) [dTZ)"^] in terms of TZ. We 
define 



5L2 
5TZ 



^ia^QTZ) - aQc^d^TZ 
at 



(28) 
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Then the curvature perturbation obeys the equation (5£2/'572.|i = at hnear level, i.e. 

^ (a^QTZ) - aQcld'^n = . (29) 

We write TZ in Fourier space, as 

7^(r, x) = — / dPkUiT, fc)e*'^-^ , 7^(T, k) u(t, A;)a(fc) + u*{t, ~k)aU~k) , (30) 
(27r)'^ J 

where a(fc) and at(fc) are the annihilation and creation operators, respectively, satisfying the commutation relations 

[a(fei),at(fc2)] =(27r)3,5(3)(fci-fc2), [a(fci), 0(^2)] = [a\ki) , a\k2)\ =0. (31) 

Note that r = J a^^dt is the conformal time, which is given by t = — l/(aiJ) in the de Sitter background. The 
asymptotic past and future correspond to r — >■ —00 and r — >■ —0, respectively. 

Let us derive the solution of TZ during inflation at linear order. The equation for the Fourier mode u is given by 

(a^QY 

u+^^^u + cl^u = 0. 32 

In the large-scale limit (k 0) the solution to this equation is u — C1+C2 J{a^Q)^^ dt, where ci and C2 are integration 
constants. Provided that the variable Q changes slowly in time, u approaches a constant after the perturbations exit 
the Hubble radius (cgk < aH). Introducing a field v = zu with z — a\J1Q the kinetic term in the second-order action 
(1221) can be rewritten as J drd'^xv''^ /2, where a prime represents a derivative with respect to r. In other words, v is 
the canonical field that should be quantized. Equation ([32]) can be written as 



cik^-—jv = Q. (33) 

In the de Sitter background with a slow variation of the quantity Q, we can approximate z" / z ~ 2/t^. In the asymp- 
totic past (fcr — >■ —00) we choose the Bunch-Davis vacuum characterized by the mode function v = e~'^'^'^'^ j \f2csk. 
Then the solution of Eq. ((33|) is given by 

--(-'^)= 2(c.fc)3/^7Q ^^ + --^-)- 

The deviation from the exact de Sitter background gives rise to a small modification to the solution (|34|) . but this 
difference appears as a next-order correction to the power spectrum and to the non-Gaussianity parameter [28| . 

The two-point correlation function, some time after the Hubble radius crossing, is given by the vacuum expec- 
tation value (0|7?.(r, fci)7^(T, A;2)|0) at r « 0. We define the power spectrum V-R,{k\), as (0|7?.(0, A;i)7?,(0, A;2)|0) = 
(27rVfci)'P7j(fci) (27r)35(3)(fci + k^). Using the solution dMl), we obtain 

" %nWs " ^^^M^FesCs ' ^^^^ 

where we have used defined in Eq. (f26| . Since the curvature perturbation soon _approaches a constant for Cgk < aH, 
we just need to evaluate the power spectrum psp at Cgk = aH during inflation j65|. The spectral index of TZ is given 

by 

dlriVn 



dink 



where 



Csk—aH 



~2e - Sf - Vs - s 
-2es -~Vs-s-8Si + 25gx , (36) 



„ , ^^77^- (37) 
He. He. 



Here we assumed that both H and Cg slowly vary, such that d In fc at fc = aH is approximated as d In fc = d In a = Hdt. 
In the second equality of Eq. (IM|) we used Eq. (P7)) to convert e to at linear order. 
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Let us also derive the spectrum of tensor perturbations generated during inflation. In addition to the three- 
dimensional metric a'^{t)e'^^5ij coming from the scalar part, we consider the intrinsic tensor perturbation h[J\ In 



terms of the two polarization tensors and ef- we can write h 



(T) 



h+elj + hxCij, where both the symmetric 



tensors Cij are transverse and traceless. We also impose the normalization condition, eij{k) eij{—k)* = 2, for each 
polarization, whereas efAk) e^A—k)* = 0. The second-order action for the tensor modes is given by 



(38) 



where A = +, x, and 



2 _ t«4 



My 



1 + 8(5^ + 0(e2). 



(39) 
(40) 



A canonical field associated with h\ corresponds to vt = zxhx and zt = a^/2QT- Following the same procedure 
as before, the solution to h\ recovering the Bunch-Davis vacuum in the asymptotic past is 



hx = 



2(cTfc)3/2VQ: 



(1 -I- icrkr) , 



(41) 



which approaches h\ — iH /[2{cTk)'^^^ y/Qr] after the Hubble radius crossing. According to the chosen normalization 
for the tensors e^^ , the two-point correlation function leads to a tensor power spectrum, whose expression is given by 
VT^^-k^\hx\y{27r^),i.e. 



Vt 



2H' 



27r2 Qt4 7r2M2F 



(42) 



where in the last approximate equality we have used Eqs. (15^ and pOj) at leading order. The spectral index of Vt is 

dln-Pr 



ut 



dink 



CT^—aH 



= -2e-(5, 



= -2e, -8(5f 



2(5, 



GX 



(43) 



which is valid up to first order of slow-variation parameters. 

For those times before the end of inflation (e ^ 1) when both Viz and Vt remain approximately constants, we can 
estimate the tensor-to-scalar ratio as 



Vn 'Qt4 



IGcce.,, 



(44) 



where in the last approximate equality the terms at second order of slow-variation parameters are not taken into 
account. Using Eqs. and (jH]), we obtain the consistency relation 

r = 8c, i^UT - 86^ + 26gx) ■ (45) 

The Gauss-Bonnet term, as well as the Galileon term, modifies the consistency relation r = —ScsUt valid in 
k-inflation iGSl. 



IV. THIRD-ORDER ACTION AND PRIMORDIAL NON-GAUSSIANITIES 



In order to evaluate the three-point correlation function of curvature perturbations, we need to expand the action Q 
up to third order in the perturbation fields. From the computational point of view, the fact that we just need to 
expand the lapse and shift only up to first order (as already stated in Sec. IIII[) is very useful for deriving the action 
at cubic order. 
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After many integrations by parts, the third-order action following from ([2]) can be written as 

S3 ^ J dt d^x {ai + (as 7e + 03 7t + 04 d'^U/a^ + ag^V/a^) 

+ a [ae d^TZd^^/a^ + 07 7^7^ + ag nd^Tl/a^ + ag [didjipdidjip - d'^ipd'^i)) / 

+ awididjipdidjU - d^-^d^V^ila" + an 7^aV/a^ + axiild^i)la^ + ai^lld^n/a^ + ai4 {dTZf/a^ + aisU^] 
+ bi + b2n{dnf/a^ + b3n^n + ci n^^n^,^/a'^ + C2n^d'^'ip/a'^ + C3nnd'^^/a^ 

+ [diil + d27^) {d.djijdidji; - 9V9V)/a^ + rfaftTe^.V^V/a^} , (46) 

where 

ai = UI^^FH^ + iMl^FH - XPx - 4:X^P.xx ~ 4:X^P.xxx /3 - SOH^i 

- 2H<p{10XG,x + nX^G,xx + 2X'^G^xxx) + 2x0,4, + WX^C^x/S + 4X3G,0XJf /3 , (47) 
02 = W3 = ~9M^iFH^ ~ 9M^iHF + 3{XP,x + 2X^P,xx) + lAAH^i 

+ 18iI0(2XG^x + X'^G,xx) ~ 6(XG,0 + X^G^^x) , (48) 

(49) 
(50) 
(51) 
(52) 
(53) 
(54) 
(55) 

In Appendix [B] we present more details for the derivation of Eq. (P^. Terms of the forms 'R? and aR? vanish because 
of the background equations of motion, which is similar to what happens for the terms 1^} and aR, in the second-order 
Lagrangian ([T^ . Another analogy between the two actions Si and 5*3 consists of the cancellation of the self-coupling 
term for the field V': that is why the cubic term in ■0 is absent in Eq. (1461) . In addition to the curvature perturbation 
72. the action pS)) involves the terms a and ■)/'■ Using the constraint equation pT| to eliminate a, the action pS)) 
reduces to 

5*3 = y" dtd^xa^iA^il^ Aik^d'^Tlla^ + A3V}d'^i;lc? -(- + {A^Tl + A^Tl) {d^dj^d.d^^ - d'^^d'^ip)/ a'^ 

+ Aj'R{d,dji}d,dj'R - d^iljd^'R)/a'^ + A^-Ri&Rf/a^ + A^djld.^l^ d^i^ / a'^} , (56) 

where 

Ai = bi + Liaiz + L\a3 + L\ai , A2 = Li(Liai + a^,) , A3 = C2 + Lia^ + Lfa^ , 

Ai = 63 + Lifly + Ljai , A5 = Liag + di , Aq ^ di , 

A7 = Liaio, As = b2 + ai3Li/2 + Li{di3 + Hai3)/2, Ag = d3 . (57) 

For the derivation of Eq. ([55| the term proportional to TZTZd^TZ has been integrated by parts, as 

c{t)nnd^n = c{t)n{dnf/2 - c{t)n{dnfl2 + total derivatives . (58) 



03 


= -Sag = -3[2M^ii^iJ M^yP - 48H^i - 20(2XG,x + 


- X'^G^xx)] , 


04 


= -IGHt 




ae 


= -a7/9 ^ an = -~W2 - ~-[M^y{2HF + F) - 2<j)XG,x - 


^ 24ij2^] , 


as 


= 2aio = 2bi = -2c2 = -4di = 16^, 




ag 


= ai2/4 = -ai5/6 = -(Af^^^^ _ 2AHi)/2 , 




ai3 


= 2ai4 — 263/9 = — ci = — C3 = —4^2/3 = ^3 = —2wi = 


-2(Af2jF - 8i7e) , 


&2 


= W4 ^M^yF-8^. 





The next step is to eliminate ip in the action ()56p by using the relation ()25p . In doing so, we introduce a new 
perturbation variable X, as X = wix/a^ (in which case d^X = QTZ). This definition is suitable in order to obtain 
the contributions proportional to the linear differential equation of TZ, coming from terms which are in the form d^X. 
Substituting the relation ip — —LiR. -\- a?X /wi into the action (l56l) . it follows that 

S3= fdt d^x (aVi + a/2 + h/a) , (59) 
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where 

A 



A, + A,^-AX]Ti' 



Wi 



+ — [AsU + AeU) {^^^JX)i^^^JX) 



2LiQ 



Wi 



Q_ 

Wi 



A2 - A3L1 + A5^^^^ ~ Aj^ ] v}d^n + Ai^:^^niid^n + A^mmf ~ AQ^^ii{dnf 



Wi 



Wi 



n{d,djn){d^djX) - d^nd^ndiX , 



Wi Wi Wi 

2 A T \ lb !fa a T>\l a a T>\ /'«2t^\2i I A t'^T) r/a a t>\I i^t t>\ /'«2t^\2i I 4 r2/ 



(60) 



(61) 



h = {A5LI - ArLi) n [{d^djn){d,djn) - [d'^Tlf] + AqLItI [{d^djTVtid.djU) - [d'^Uf] + AgLUdUfd^TZ .(62) 
In standard inflation with a canonical kinetic term (i.e. F = 1, P = X — V{<j)), ^ = 0, G = 0) one has Ai = —M^^e/H, 



MlJ{2H), As = 3Af2^/2, A 



A2 = A3 = Ay - 0, A4 = SMpV, As = 

and wi = M^j, where e = (j)'^ / {2H^ M^^) . It then follows that 



pi/ 



-MpV, Ag = -2M2j, ii = 1/H, Q - MpV, 



/i = Af^Je(l-e/2)7^2(37^-7^/i^)-2e7^a,7^(9,A' + (37^-7^/i^)(9,aJA')(9,aJA')/( (63) 
/2 = -M2le(7^ - 2n/H){dnf + {M^ie/H){3n - il/H)nd'^H - (l/i^)(37^ - n/H){d,djn){d,djX) 

+ {2/H)d^nd,nd^X , (64) 

/3 = MpV(2i^2)(37^ - 7^/i^) [(a,5J7^)(9,a,7^) - (927^)2] _ {2Mly/H^){dnfd^n , (65) 

which agree with Eq. (114) in Ref. [131 in units where Mpj = 1. 

Although the expressions for the building blocks of the third-order action are correctly given in Eqs. (|60p - (p^ . it is 
more convenient to perform several integrations by parts to bring each contribution into a simpler and more usable 
form. In Appendix [C] we show that the integrand /3/a in Eq. ([CTl) reduces to 



/3/a = (qa/a) {{dTZy (d^TZ) - TZ d,dj [{d{R) {djTZ)]} + total derivatives , 



where 



03 = Afiif 



2 ad/ A^L\ — AyLi 



(66) 



(67) 



Along the same lines the term a/2 can be written as (see Appendix [D] for details) 

aq2 



a/2 = - 



where 



■(a^g7^) {(a7^)^ - d-'d,dj[{d,n){d3'^)]} 

2aLiQnnd^n + a (A2 - A3L1) U^d^U + a 



As + - — iaLiQ) 
a at 



{{d^x) {d,n) {d'n) -nd,dj[{diTZ) (d^x)]} 

n{dnf + total derivatives , (68) 



q2 



AAqLi 2 d / Aj - 2A5L1 
wi dt \ a^wi 



2A9L1 



Finally, the term a^/i, as shown in Appendix lEl is equivalent to 

d 



A4 + qi{Q + 'SHQ) - QgiJ 7^7^2 - 2ql7^7^ —{a^QTZ) 



A1+A3— -glg)7^^ 

Wi 



-,3 r 



Wi 



a' 








3HA5' 


2 


wf 


--( 

dt ^ 




+ 2^ 



. Q , QAe „ d 

Ag 1 Qwi — 

wi wi at 

2A5 d 



Wi 




(4)' 


1 -SQHA,} 




Wi 



(69) 



ndiUd^X 



{d^'R){dxf - ^(«'Q7t) {{dkn){dkX) - d-^d,dj[id,n){d,x)]} 



total derivatives . 



(70) 
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In what follows we omit to write the terms corresponding to total derivatives. We determine the coefficient qi such 
that the term -2qiTZ'R. -^{a^QTZ) in Eq. is merged with the term -2aLiQmZd^TZ in Eq. dBH]) to give rise to a 
combination proportional to the linear equation of motion for TZ [i.e. Eq. (f29| ]. This can be achieved by demanding 



-Li/cl 



Since, in general, we have that 



c{t) {{d,x) {d,n) (d^TZ) -nd,d,[{d,TZ) {djX)]} = c{t){{d,x) {d,n) {d^n) ~ {d-^d^n) d,dj[{d,n) {d,x)]} 

the last term in Eq. ((70)) can be merged with the term including <?2 in Eq. (|68)) as follows: 

A. 



(71) 



(72) 



{{dkn){dkX) - d-^d,dj[{d,TZ){d,x)]} 



^^{{dkTZ){dkX) - d~'d,d,mnMX)]}^-^ 



92 _ 2c|A5Q 
2 



id^n){{dkn)idkX) - d-^d,d,mn)id,x)]}, 



(73) 



where S£2/STZ\i is defined in Eq. 



Along the same lines one can show that 



c{t){{dTZf {d^TZ) -nd,dj[{d,n) [d^n)]} = c{t){{dnf (d^n) - (d-^d^n) d^dj[{d^TZ) {djTZ)] 

= c{t)d^n{{dnf - ^-^^^^J[{^^n) (djlZ)]} . 



(74) 



We use this result to merge two different contributions, namely the term /3/a in Eq. (|66p and the first term on the 
r.h.s. of Eq. §8^. We have 



1 



--^ (A7 - 2^5X1) ^U^Qil) + aq^d'^n 
Zwi at 

1 



1 



sn 



93 - ^ (A7 - 2A5L1) 



Finally, merging all the contributions together, the third-order action (|56p can be written as 

where 
C3 - 



(75) 



(76) 



d^a;|a•^ClMJ^l7^7^^ + aC2M^{R{d'R:Y + a^C3Mpl7^•^ + a^C47^(a^7^)(9,A') + a\C5/M^^)d^n{dXY 

aCQ-R^d^n + (Cy/a) [d'^n{dnf - nd^dj{din)(djTZ)] + a{Cs/Mpi) [d'^TZdMX - ndidj{din)(djX] 
SC2 



5n 



(77) 
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The dimensionless coefScients Ci {i = 1 
1 



are given by 



Co = 



1 
1 



A4 + qi{Q + 3HQ)-Qqi 

1 



_Q_ 



3- 



_Q_ 

HQ 



dt 



Li 



a at 

Q 



A/pi 
Q_ 

Wl 



= Ml + ^3— -qiQ] = —<L 



Ml, 
1 



1 d 



Ml,F-8(+--{aL^{Q-Wi)) 



Wl 



— {Ae + Ag)~wi — 
Wl at 



pi 



Li{Liai + as) + 015 + (ai2 + Lias)— + % 

Wl Cg 



( 1 - - 



A. 



3HA5 



Wl 



_Q_ 

2wi 



1 + 2wi 



A, 



3H- 



.A. 



,t2 
Pl 



Aa 



d (A 



WX 



dt \ w\ 



mA^ 



wX 



Ml 
2w\ 



-M'F{1 - HLi) - 2AHS, 1 - -HLi 



2dt \wl 



pl : 



Li{2MlyF - Lias - 64^0 + 24^ 



Ce = A2 - A3L1 = Li 

C7 = 93 - ^ {A7 - 2A5L1) - III \mIiF{1 - HLi) - SHiiA - 3HLi] 
2wi 6 L ^ 



2wi' 



Li{MlyF -24.Hi) + lQi 



AL 



pl 



92 _ 2clAzQ 
2 w? 



+ (l? [Li(A/p2 L^ - 24i70 + 24^ 
^(Afp2 - 24H0{HLi - 1) + ^ 

Wl ^ Wf 



(ft 



Li[Li(Af2F-24LrO + 16e] 



2wi 



Li{M?^F-24H0+8i 



Mpi , 



(78) 
(79) 
(80) 
(81) 

(82) 
(83) 

(84) 
(85) 



with A5 = -{Li/2){Ml^F - 2AHi) - 4^. Note that we have used A^ = 3Q to simphfy Ci. The coefficient in front of 
5C2/5'R\i in Eq. dTT]) is 



Ti = ^ {{dk'R)(dkX) - d-^d^d,[{d,n){d,x)]} + qiun 



A7 - 2A5L1 

Awio^ 



{{dnf - d-^d,d,[{d,n){d,n)]} . (86) 



Each contribution to J-i includes spatial and time derivatives of 7?., which vanish in the large-scale limit {k — > 0). 
Moreover the term 5C2/5TZ\i survives only at second order in TZ. When we evaluate the three-point correlation 
function of TZ below, we neglect the contribution of the last term in Eq. ([77]) relative to those coming from other 
terms. In Refs. [l^, [13, HI] the term proportional to is present in the definition of Fi. After a suitable field 
redefinition this gives rise to a correction of the order of the slow-roll parameter 77 — e/{He) in standard inflation [6^. 
We have absorbed such a contribution to other terms in Eq. ([77| [such as Ci] and hence the field redefinition is not 
required in our method. 

The Hamiltonian in the interacting picture is given by Hint = —^3 [SH^jll^- The vacuum expectation value of 
TZ for the three-point operator at the conformal time t — Tf can be expressed as 



(7^(fcl)7^(fc2)7^(fc3)) 



dra (o| [niTf,ki)niTf,k2)n{Tf,k3),nintiT)] |o) , 



(87) 



where is the initial time when the perturbations are deep inside the Hubble radius. Since r ~ — l/(aiJ) during 
inflation it is a good approximation to take —00 and Tf 0, where the latter corresponds to some time after 

the Hubble radius crossing. 

In order to evaluate the vacuum expectation value (j87p we use the curvature perturbation (j30l) in Fourier space 
with the mode function u{T,k) given in Eq. (1341) . Each term in the third-order Lagrangian (j77p includes the phase 
factor of the form / ^^3^ gi(fc4+fc5+fc6) :c^ which gives rise to the delta function (27r)^ (5^'^^ (fc4 + + fcg)- Among the 
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combination of the annihilation and creation operators, the non-vanishing components are 

{0\aiki)a{k2)a{k3)aH-ki)a\-k5)a^{-ke)\0) = {0\a{ki)a{k5)a{ke)a\-ki)aH-k2)a\-k3)\0) 
= (27r)^((5(3)(fc4 + fci) 



S^^\k5 + k2) S^^^ke + fcs) + S^^^\k5 + fcg) 5(=^^(fe6 + fca) 

+ (5(3) (fc4 + fca) k^^fcs + fci) 6'^^Hke + k^) + S^^^k^ + fca) S^'^^ke + ki) 
+ S^^^ki + fcg) M^^fcs + fei) ^^^^fce + fca) + (5(3) (fes + fca) (J^^Hfee + fci) 



When we carry out the integral in terms of r we assume that all the terms Ci {i — 1, • • ■ ,8) slowly vary in time 
relative to the scale factor a, so that it is a good approximation to treat them as constants for the integration. If this 
condition is not satisfied, one needs to numerically solve the integrals without approximations. In the following we 
present each contribution of the three-point correlation function coming from the integrals given in Eq. (j77p . 



• (1) H^:^ = 



J d^xa^CiM^{R'R: 



{n{k,)n{k2)n{k:,))^^^ = {2TTf 5^^\k, + fc^ + ^3)^^ 



1 



K 



if2 



sym 



16Q3c6 {kik^hf 

where K = ki + k2 + k'i. The symbol "sym" means the symmetric terms with respect to ki,k2, k^ 
. {2)H^,^i = -Jd^xaC2M^,n{dnr 

C2Ml,H^ 1 



(7^(fcl)7^(fe2)7^(fe3)) 



(2) 



(fci • fca + fcs • fca + fea • ki) ~K + 



16Q3c8 {kik2hY 

kik2 -f fc2fc3 + k^ki kik2kz 



i3)H, 



(3) 
int 



-/ d^xa^CaMpiTZ^ 

(7^(fel)7^(fe2)7^(fc3))(3) = {27rf 5^^\ki +k2 + kj] 



3C3Afpiff5 



1 1 



(4) Hi 



(4) 



(7^(fel)7^(fc2)7^(fc3))(^' = (2^)35(3)(fci + k2 + k^) 



32Q2c6 {kik2k^f 



(fcl • fc2)fc3 

K 



fcl + fc2 



sym 



• (5) ^-Jd^xa\C,/Ml)d^n{dXf 

(7^(fcl)7^(fc2)7^(fe3))(5) = (27r)3,5(3)(fci + fe2 + fc3)- 



16QMp2c6 (fcifc2fc3) 



kl{k2-k^) 



K 



1 + 



sym 



(6)77. 



(6) 
int 



Jd^xaCelZ'^d^n 



(7) 



(7^(fcl)7^(fc2)7^(fc3))^^^ = (27r)35(3)(fc^ + fc2 + fca 



/d^x (Cy/a) [a27^(a7^)2 - 7^a,9J(a,7^)(9J7^)] 



1 1 



4g3cf fcifc2fc3 



(7^(fcl)7^(fc2)7^(fc3))*'^ = (2^)^5(3) (fci + fc2 + ^3) 



(89) 



(90) 



(91) 



(92) 



(93) 



(94) 



1 



1 



.^3^10 (fcifc2fc3)3 if 

[fci(fe2 • fes) - (fcl • fe2)(fei • ks) + sym] . 



fclfc2 + fc2^3 + k^ki 3fclfc2fc3 



(95) 
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(8) H^^i = -jd?xa(C8/Mpi) [d^ndiUd^X - lld,d,[d{ll)[d,X)\ 



{TZ{k^)n{k2)'R{k:,))^^^ = (27r)3j(3)(fei+fc2 + fc3) 

2 



1 1 



32Q2Mpicf (fcifczfcs)^ if 

2fci + fc2 + fc3 2fci(fc2 + fc3) 



is: 



[kl{k2 ■ fca) - (fei • fe2)(fci • fca)] + sym \ . (96) 



We write the three-point correlation function of the curvature perturbation, as 

{n{ki)n{k2)n{k3)) = (2^)3,5(3) (fci + ^2 + k3)Bn{kiM, h) , 

where 

B-fi{ki,k2,k3) = ^-^^T^—^An{ki,k2,k3) . 



(97) 



(98) 



Recall that the power spectrum V-ji is given by Eq. ([35]). Collecting all the terms in Eqs. (|89l) - ([96|) . it follows that 



A 



n 



Ml\l I 2 

Q 



3 _H_ (fcifc2fc3)2 , i ^ V fc^ - i V - — V fc^ fc A 

fe) ]^ Efc^^ + ^Efc^^l-^E^^^^I-fci^^fcaE^'^: 



II Q \ 1 

4 

1 f H Y I 



Z ( H V (fcifc2fc3)2 



^'^c^Upi 



if3 



2ct VMpi/ if \ ii:^ 



3fcifc2fc3 



1 ii Q 1 
8^ Mpi if^ 



2 

l>j 



klk2k3 E - 2 ^1^2fc3if ' - 6 E _ ^ fcf + - if E fc,^ 

i i^j i i 



(99) 



which has an explicit dependence on the wave numbers. The bispectrum (j99p is the central result of our work. 
In general the non-linear parameter /nl associated with non-Gaussianities of the curvature perturbation is defined 

as 



/i 



NL 



10 An 



(100) 



3 Eti kf 

which matches with the notation of the WMAP group. For the equilateral configuration with ki = k2 = kj, one has 

2 



pcquil 



40Mp^i 



"^NL 9 Q 



Lc 

^ 96cr' 72 Mpi ' 24 M^^ ^ 24 1^ Ai^^ ^ 



1 / ii 



36c2 VAipi 



13 / ii 



17 H Q 
192^Mp^M|; 



(101) 



96c4 VAipi 

which is independent of the wave numbers. Recall that the coefhcients Ci (« = 1, • • • , 8) can be evaluated by Eqs. ([78 



V. EXPANSION IN TERMS OF SLOW VARIATION PARAMETERS 

Let us derive a simpler form of /nl"'' by using the approximation that the slow- variation terms defined in Eq. ([5]) 
are much smaller than 1. We also use the expansion of Li and eg given in Eqs. (I^IH) and ([?7)) . respectively. First we 
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write the variables Q and in the forms 

4Mp[E 2 2w'lw2H — w|w4 + 'iwiWiW2 — 2w\w2 



where 



r. "-^— . (103) 



We also define the following quantities 

X = F^ [x^Pxx + 2X''Pxxx/i + iH{XG,x + SX^G.xx + 2X''G^xxx) - 2{2X^G^^x + X^G^^xx)li\ , (104) 
and 

Xg^XG^xxIG^x- (105) 

It is convenient to introduce S and A in order to obtain a cornpact expression of /nl"''- Note that these parameters 
have been already introduced in the context of k-inflation [27l . [28| . 

We first replace with respect to E, which in turn can be written in terms of e^. Using Eq. (j26p . the parameter 
ts can be expressed in terms of other slow- variation parameters. Then we replace F as F = QcH (e^Afpi)- Finally we 
expand the coefficients Ci {i — 1, - ■ ■ , 8) in terms of the slow- variation parameters, by keeping and \g as unknown 
parameters. A similar method is applied to all the other terms, except for C3. In this case the terms Px and G^xxx 
are expressed in terms of E and A, respectively. We replace any left G_xx with Aq. We also multiply and divide A by 
E, and we consider the ratio A/E as a free parameter of the theories (as we do for and Ag). Then, as done for the 
other terms, we replace E in terms of Q, and any left F in terms of Q and e^. Finally we use the relation between e 
and other slow- variation parameters. By following this procedure, the final expression of each /^l^'^'*^ (i = 1, • • • ,8) 
only depends on e^, r]s, s, Sp, S^, Sgx, and the free parameters c^, Ac and A/E. This form allows an expansion with 

respect to the si 
Ci are given by 

/Z" - y (1 - ^) - ^ (8 + + 4 - - (106) 
/X"^^' - ^f4-lVT7^(^^ + 8^«-2. + ,.), (107) 



respect to the slow- variation parameters. Then the leading contributions to each /^l^^^^*^ coming from the coefficients 

/.cquil (1) 



108 \c4 J lOScj 

8iU-0-8Te+T62^(^'^ + 2^---^-)-16-2 
_5_ Sgx 
27 

10 

27 ct 



/Nf'^'^ = i,{^-lj-f^^ + j^^{SS, + 25Gx-5,)-^{6SGx^S, + i\GSGx+SS,) 

ci{l + XG){6F~SS^~2SGx) , (108) 

^e,uU(4) ^10^^ ^^^g^ 



^cquii(6) _ 20 (1 + Ag) Sgx nm 
^^i^ " 81 7s ' ^ ' 



nequii(7) _ 65 Sgx , „^ 



.cquil (8) _ 85 Sgx . „^ 

/nl - ^2 ■ ill«^i 

Let us now add up all the contributions together. In doing so we take the largest contributions, that is, we discard 
the corrections of the order of slow-variation parameters relative to other existent terms. For example, if a term 
Sgx / i^scl) already exists, we ignore the terms like Sgx/c^- Then we have that 



pequil 85 1 

/nl - 324 I ^ 





10 


A 


55 es 


5 


85 s 


c? ) 


81 


E 


36 c? 


12 c? 


54 c? 
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where we have tacitly assumed that < 0(1). The result (|114[) is valid for any quasi de Sitter background. 

VI. CONCRETE MODELS 

In this section we apply the results derived in previous sections to concrete models of inflation. This includes (i) 
k-inflation, (ii) generalized Galileon model, and (iii) models based on BD theories. 

A. k-inflation 



Let us first consider k-infiation models with F = 1,^ = 0, and G = 0. Since Cg = e = —H/H^ = XP^x /{M^iH^) and 

Li = 1/H in those models, it is possible to derive an exact expression of /^l"'' without employing the approximation 
in terms of slow- variation parameters. The coefficients (|78p - (j85p are given by 

C3 = - 



Ci = 4(e - 3 + 3c2 - ry) , C2 = ^{l - + e + - 2s) 
dt ci 



pi 



2A 



C4 



C5 



1 



Cfi — C7 



C8 



0. 



(115) 



where r] = r/s ^ ^/{He), S = XP^x + 2X^P,xx and A = X^P 



XX 



2X3Pxxx/3. 



The field propagation speed squared is 



P 



.X 



XP 



X 



(116) 



fcquil _ ^5 ( -1 ^ 1 _ 

/nl ~ 324 I ^2 I 



10 A 55 e 5 ?7 85 s 
81 S ^ 36^1 12 Z ~ 54 c! ' 



(117) 



" P.x + 2XP.XX S 
Plugging the coefficients (jll5l) into Eq. (|101l) . it follows that 

1 

where we have ignored the second-order term /{2c^^) in the expression of C4. The non-linear parameter (jll4p 
derived under the slow-variation approximation also gives the same result. Note that this coincides with the result in 
Refs. [13, HI. 

In standard inflation driven by a potential energy V{(()) of the field 0, i.e. P — X — ¥{(()), one has = 1, A = 0, 
and s — Q. Then the equilateral non-linear parameter (jll7p reduces to 



/•cquil 

/nl 



55 
36* 



12 



(118) 



which means that 1/^1"'' I ^ 1 



Generalized Galileon model 



equil 



Let US consider the generalized Galileon model with F — 1 and ^ = 0. Since the exact formula of f-^^ 
complicated, we employ the result (|114[) derived under the slow-variation approximation: 



IS 



- 324 

Using the background equation Q , the term Sqx / can be expressed as 

^GX ^GX H4>G,x 



85 





10 


A 


55 eg 


5 7?^ 


85 s ^ 5gx 


■20(H-Ag) 65 




81 




^ 36 c2 ^ 


12 d 


54 cl es 


81 ' 162c2 



e + 5gx Px + (4i/(/) - (j))G^x - 2G,, 



In the limit that cj^l the dominant contribution to /^l"'' is 



j^equil 
'NL 



85 1 10 A 



324 c? 



81 E Px + (4H(t) - <j>)G^x - 2G,, 



20(1 + Ag) 
81 



65 



162c2 



(119) 



(120) 



(121) 
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This matches with the result of Ref. in which the authors ignored the terms 4>G^x and G^^ relative to H^iG^x- 
The last term in the parenthesis of Eq. (|12ip gives rise to an additional contribution to the terms appearing in 
k-inflation. 



C. Brans-Dicke theories 



The action in BD theory with a field potential V{(j)) is given by 



S — I d'^x\/—g 



(122) 



where wbd is the BD parameter. Compared to the original paper 6l[ we have introduced the reduced Planck mass 



Afpi, so that the field </> has a dimension of mass. In this theory the background equation ([7]) gives 

2ei = -62 + e2e3 + WBD eL (123) 

where 



H 



ei = --2, £2 = ^, £3 = ^. (124) 



Let us derive a full expression of /^l"'' without using the slow- variation approximation. From Eqs. (j23|) and (|24p it 
follows that 

_ (3 + 2a.BD)£i^Mpi 2_ 

The time-derivatives appearing in the coefficients Ci can be expressed in terms of (i = 1, 2, 3) by using the relation 

£2 = ffe2(ei-e2 + e3). (126) 
Then the coefficients Ci (i = 1, • • • ,8) are given by 

(3 + 2a;BD)el [(2a.BD-3)e2+4(e3-3)e2 + 8e3] 
(3 + 2c^BD)ei[8e3 + (4e3-12)e2 + (3 + 6c^BD)ei] 

^^^^ i^TW ' ^ 

C3 = Ce = Ct = Cs = , (129) 
(3 + 2ujBD)el [-16 - 16e2 + (2c^bd ^ l)^!] 

^ ^(^T^^)^ ' ^'^^^ 

Afpi (3 + 2u;BD)el 
4(2 + e2)2 ' 



which lead to 

^equii _ 5 4863 - 72(1 - e3)£2 - (6 - 68c^BD - 36e3)ei + (48 + 68c^bd + 6e3)e^ + (12 + II^bd - 2c^iD)e| 

"18 (2 + 62)4 • ^ 

Since |e2| ^ 1 and |e3| ^ 1 during inflation, we have 

/^r'^-^62 + ;^63, (133) 



/■oquil 5 5 



which gives I/^l"''! ^1- It is worth mentioning that the same result as Eq. (jl33l) follows by using the approximated 
expression of /^'l"''''' (* = 1; ' ■ ' ^8) derived in Sec.|Vl The above analysis covers the case of f{R) gravity in which 
the BD parameter is ojbd — 0. Hence, the Starobinsky's inflation model, f{R) = R + /{6M^), leads to 1/^1"''! ^ 1 
as in standard inflation. 
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D. Potential driven inflation in the presence of the nonminimal couphng and the Gauss-Bonnet term 

The result (jll4p shows that, in the hmit ^ 1, the terms coming from the nonminimal coupling and the Gauss- 
Bonnet coupling are suppressed relative to the term proportional to 1/c^. The effect of those couplings on /^l"'' 
appears indirectly through the change of c^. 

Let us consider the following action 



d x\J —g 



Ml 



(134) 



which corresponds to inflation driven by the field potential V{(j)). Expanding the scalar propagation speed (j24[) up to 
second order, we obtain 



2_ 2S^{5F-Sd^){3SF~24:6^~46px) 

OPX 



(135) 



where Spx — ^^/{MIH'^F) is another slow- variation parameter. Since \Sf\, \^px\} ^ 1, the scalar propagation 
speed is very close to 1, i.e. \cj. — 1| = C'(e^)- From Eq. (I114p we have 



/.cquil 55 



12 



10 , 



3?, 



(136) 



and hence the non-Gaussianity is small. 

The non-Gaussianity can be large for the kinetically driven inflation in the presence of the terms P{4>,X), £,{4>)G, 
and G{(f>, X)D(j). We will study such models in a separate publication (68| . 



VII. CONCLUSIONS 



For the general single-fleld models described by the action ([2]) we have calculated the three-point correlation function 
of curvature perturbations TZ generated during inflation. This covers the inflationary models motivated by low-energy 
effective string theory, scalar-tensor theories, and Galileon-inspired gravity. The Gauss-Bonnet term £,{4i)G and the 
generalized Galileon term G{(j), X)[2(j) give rise to the scalar propagation speed Cs different from 1, as it happens 
in k-infiation. These models can lead to large non-Gaussianities of primordial perturbations detectable in future 
observations. 

Using the ADM metric in uniform-field gauge [Scj) = 0), we have expanded the action ([2]) up to third order in the 
perturbations. The second-order action associated with linear curvature perturbations is given by Eq. (|22[) . where 
both Q and cj. need to be positive to avoid the appearance of ghosts and Laplacian instabilities respectively. In 
the quasi de Sitter background we have solved the equation for curvature perturbations and derived the scalar power 
spectrum and its spectral index. The tensor spectrum as well as the tensor-to-scalar ratio is also evaluated to confront 
the models with observations. 

After the derivation of the third-order perturbed action ()46p . we have expressed it in a more convenient form for 
the calculation of non-Gaussianities by making a lot of integrations by parts. The coefficient J^i in front of 5C2/5TV\i 
includes only the derivative terms of 72., which vanish in the large-scale limit. Moreover the term 5C2/5TV\i survives 
only at second order of perturbations. When the vacuum expectation value of TZ for the three-point operator is 
evaluated, we just need to compute the contributions coming from the 8 terms in Eq. (j77p by treating the coefficients 
Ci (i = 1, • • • , 8) as slowly varying parameters relative to the scale factor. Finally we have derived the three-point 
correlation function as in the form (j99|) . For the equilateral configuration (fci — k2 = fca) the non-linear parameter 
/nl"' reduces to Eq. (|101l) . which does not have any momentum dependence. 

Under the approximation in terms of slow-variation parameters we showed that /^l"'' reduces to a simple form 
(I114p . Provided that <C 1, one can realize large non-Gaussianities with I/^l"''! ^ 1- For ^ 1 the terms 
involving Sp and i5{ are sub-dominant contributions with respect to the first two terms in Eq. (|114l) . but the last 
term in Eq. (jll4l) . coming from the Galileon-type self interaction, can be comparable to the dominant contributions. 
Compared to k-inflation the presence of the terms S,{(t>)G and G((/), X)[I\(I) leads to different values of Cs, so that /^l"' 
is subject to change. 

We have applied our results of the equilateral non-linear parameter to a number of concrete models of inflation. 
In k-inflation and in the generalized Galileon model, our formula of /^l"'' reproduces the previous results known in 
literature. In the inflationary models based on Brans-Dicke theories, including f{R) gravity, the non-linear parameter 
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is of the order of slow-roll parameters. This is associated with the fact that the scalar propagation speed is unity in 
those models. 

It will be of interest to distinguish between various inflationary models observationally by using our formula of 
/nl"' ^^^1 ^^'^ ^- Especially the potential detectability of non-Gaussianities by the PLANCK satellite may 

open up a new opportunity to approach the origin of inflation. 
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Appendix A: Explicit expression of Q and 

Using the expressions for wi, W2, W3 and W4, it is possible to write down explicitly the form of Q defined in Eq. (|23p . 
Namely we have 

Q = ^ [3 Mpi^F^ - 2 ^HWfG,c^x - 4 Mpi^FG,4,^^ + 2 P^x^^Mpi^F - 48 ^''H^G,xxi 
W2 

+ 2 ^^P^xxMpi^F + 12 4)^G^xMpi^FH + 6 ^^HG,xxMpi^F - 48 Mpi^FH^^ 

- 6 Mpi^F^^G.x - 48 H'^i^^G^x + 16 ^'^H^G,4,x - 16 P^x^'^H^ 

- 16 P,x4y'Hi + 192 _^ 32 cj)^HiG^^ + 3 ^^G,^^ - 16 <i>^P,xxHi] . (Al) 
Along the same lines the speed of propagation given in Eq. ([M]) can be written as 



^2 

where 

L2 = {M^iF-8Hi^^)[3Mpi^F^^^ -2AWF4>^G^^x-4.Mpi''FG,^ + 2Mpi''FP,x 
+ 2 Mpi^F <j>^P,xx + 12 Mpi^F 4,HG^x+& Mpi^F(l)^HG,xx - 48 AWF,^H^(,^^ 

- 6 Mpx^F^^^'^G^x - 48 H^i4.4>^G,x + 16 <i)^Hi^G^^x - 16 Px^H + 192 H^^^ 

+ 2,2^H ^,^G,^ + 3 </.4G,x' - 16 d^"" P,xxH i,^, - 48 ^^H^G^xx^] . (A3) 
L3 = 1536iJ4^/ + 64Oij2^/^2g,x + 128i/2^,/04G,xx + 2Mpi^F2(/,2G^^ 

- 2,MMpx^i/F^^H^ - A8Mpi^i^^''F^G,x + 24 (/.^G^x'C.^ + ^Mpi''F^G^x 

+ 24 A'WF^^i^ - 32 Mpi^F H i^^^G^xx - 64 Mpi^F H ^,00 G^x , (A4) 
Li = 8 G^x'^Uc^ 0' - 48 Gy-Ht,^ 0^ + [32 i^G,xG,^ - aWfG^ + 256 UG,xH^U4> 

- IQMpx^F^^G.xU* - l(!>UG,xP,x + l28H^S„^''G,^x ~ 16 Mpi^^^^G.^F^^] ^'^ 

- 32i^H {Si^G^xH^ + Mpi^FG^^x -iMpx^F^^G^ (/.^ 

+ [-ll2Mpx^FH^^^^G,x + 16Mpi2$^0i^,0Pjf +8Mpx^F/i^^^ - 2 Mpi'^FF^G^x 

- 256Mpi^i^F^H^i^^ + 16Mpi*i^F^F^^^ ~ 128 Mpi^ F^^ + 2Mpi*F^G,^x 

- 32Mpi^i^F^G^^ + 1536HX^^I^^] 0^ 

+ [8H (-384 i/H* + 96 Mpi^^^^F^H^ + 8 aWf^^G^^ - 6 Mpx^F/t.^ 
+ Afpi^F^G,^ - 4Mpi2FPxe0)] ^ 

+ Afpi^P (2 AWFPx - 4 Mpi^FG,^ + 192 i/^^,^^ ^ 3 Mpi^F^^ _ 48 AWF^H^^^^) . (A5) 



Appendix B: Third-order action 



In the Appendices iBl [Cl [Dl andlElwe shall use the symbol = for the quantities which are valid up to total derivatives. 
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We first perturb the action S = J d^x ^J—g L given in Eq. ([2]) up to third order. Since we choose the gauge in which 
the scalar field is unperturbed ((5(/i ~ 0), the functions P and G are expanded due to their X dependence. After the 
Taylor-expansion of the action up to third order, we can perform the following steps in order to simplify the result. 

1. We start by removing the cubic term in -0. One can employ the following relations 

c(t,x)dl^ ^ -(dic){d1^) for i==l,2,3, (Bl) 
c{t,x){dvMdod1^) = -(a;c)(%^)(%V-)/2 for i^j, (B2) 

where c depends on other perturbation variables, and a hatted index is not summed. The '(/'-cubic term auto- 
matically cancels out after the integration by parts. 

2. Now we simplify the term cubic in a. This can be done by employing the following relations 

c{t)o?a = -caV3 , c(i, x) d^a = -{di(x)(dic) . (B3) 

3. Next, we simplify the cubic term in TZ. In this case we can make use of the following relations 

c(^)7^^7^ = -c7tV3, c{t)n^1i = -6v? ji-^ciiv? , c\t,x)ii1i = -cv? 12 , c{t)v?ii = -cv? . (B4) 

Other useful relations are 

c{t)n{din)'^ = -7^ [c {diU)^ + 2ca^7^a^7^] , c{t)n^din = ~2cii{d{ii){dtii) , 
c(t)nndfn = cnidiUf + 1 [ii{dinf + 7^(a^7^)(a^7^)] + 3c7^((9^7^)(a^7^) , 

c{t, x) dfn ^ -(a^c)(a^7^) , c{t, x){din){dilZ) = -c (a^7^) 2/2 . (B5) 

After performing these integrations by parts, we find that the 7?.-cubic term can be written in the form ci {t)Ti? + 
C2{t)'Ril^ + c:i{t)'R? + C4(^)7^ {d{Rf . 

4. Now let us simplify the term quadratic in tp and linear in TZ. First we integrate by parts any derivative for 
the field 7?., so that these terms can be written as 7?. x (quadratic term in iji). Afterwards, we use the following 
relations (where i ^ j) 

c{t, x){diij)id?i:) = -c(a?V')(a|v) - c(a?v)(a|v) , c(t, a;)(%^)(%V') = -c (%V')(%V')/2 , 

c(t, x)diij - -(a,-c)(a?V) , cit, x){d,~d?i;) = ~{dU){djc) . (B6) 

5. Along the same lines we can simplify the term quadratic in -0 and linear in a after integrating by parts any 
derivative of a. 

6. Let us next consider the term quadratic in a and linear in TZ. In this case one can first integrate by parts all 
the second derivatives for the field a. Then one can integrate the following first derivative terms 

c(t, x)aa = -cc? 12 , c(t, x)adia ^ . (B7) 

7. Now we consider the term quadratic in TZ. First of all we can remove any derivative from the subset of terms 
which possess the field a. Then we eliminate the second derivatives for the field TZ by using 

c(t,x)iz1Z ^ -ciZ^ 12 c{t,x)TZiZ^ -ciZ^ + 'cn^ /2, c{t,x)iZ = -cTZ, 

c(^)(a^v-)(a^7^)(a?7^) - -c (a?V')(a^7^)V2 , c(^)(a^,-V')(a^7^)(aJ7^) - -c7^(a^J7^)(a^JV') + cTz^dfd?ij)/2 , 
c(^)7^(a?7^)(a|0') = -c{din)Hd?^) + cTz^{dld]i,)i2 , c(^)(a^J7^)(aJ7^)(a^^) = -c [dfzf{d1i,)i2 , (bs) 

where i ^ ]■ The remaining terms can be simplified by using the following expressions 

c{t,x)iZ{diiZ) = -~iZ^{dic)/2, c{t,x)n{d{R) = -7^2(a^c)/2, c{t,x)niZ = -c7^V2, 
c{t)iZ{diTZ){dii}) ^ c7^2a?V;/2 + cTZ^dji}/2 ~ cJZ{diiZ){dii}) , c(i, x){diiz) = -iz{dic) . (B9) 

8. For the terms which consist of a combination of three different fields, we can integrate by parts any derivative 
of a. 
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Appendix C: Integration of the term /3/a 



We rewrite the integrand /3/a in Eq. (j59p in a more convenient form. We define 

X = c{t)n[{d'^TZf - {didjn){d,djn)] , 



Xk = c{t) n [dkU [d^n) - ^^n{^^^kn)] - c{t) n [{dkii) {d^n) - {d,n) {d.dkn)] . 

Then, for any function c{t) dependent on t, we find 

X + 2dkXk - c7^ [{d^TZ)^ - {d,djTZ){didjTZ)] + 3c 71 [{d'^nf - {d^djn){d,djTZ)] , 
which imphes that 

c{t)Tl[{d^TZ)^ - {d,djTZ){didjTZ)] = -cTZ [(d^TZf - {didjn){d,djTZ)]/3 . 

Let us also define 

Y = c{t) {dnf {d^n) , 

Yk = c{t) [27^ djTZdkdjU - 2ndknd^n - dkU [dUf] . 

It then foUows that 



so that 



ay + dkYk = -2c 7^ [{d'^nf - {d,djn){didjn)] 



c{t) {dnf (d^n) = -2cn [{d^nf - {d,djn){d^djn)]/3 . 



By using Eqs. jCH) and JCSl, one finds 



a a 
. 1 



{A5LI - A7L1) n + AeLinl [{d,d.jn){didjn) - [d'^nf] 



A9LI 



{dnyd'n 



Employing the following relation 

c{t) TZ [{d^djTZ){d^djTZ) - [d^Tlf] = -c (d^TZ) {dj'R){didj'R) - c (diTZ) TZ {d^d^TZ) ~ cTZ [d^TZf 

= -cnd,dj[{dr1Z) {djTZ)] + c{dnf {d^TZ) , 



(CI) 
(C2) 



(C3) 
(C4) 



(C5) 
(C6) 



(C7) 
(C8) 



TZ [{d,djTZ){didjTZ) - {d^TZ)^] . (C9) 



(CIO) 



we obtain the form of /3/a given in Eq. (|66 



Appendix D: Integration of the term a/2 



Regarding the integrands of the term a/2 in Eq. (|59p , we start by considering the contribution written in the form 
c{t)iZ{didjTZ){d^djX). We define 

Z = c{t)n{didjn){d,djX) , (Dl) 

Zk = c{t){izd,n~ndrn){dkd,x), (D2) 

so that 

Z + dkZk = cn{d^djn){d,djX) + 2ciZ{d,djTZ){didjX) + cTZ {d,djTZ){didjX) + cQ{iZdiTZ - TZdJZ) d.TZ . (D3) 
This gives 

c{t) iz{^^^J^Z){^^^JX) = -[(c + ags) TZ {didjTZ){didjX) + cTZ {d^djn){didjX) + cQ (7^9^7^ - TZdiTZ) di'}Z]/2 

+ aq2 TZ {d.djTZ) {d,djX)l2 . (D4) 
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For later convenience wc added and subtracted the quantity aq2'R-{didjTZ){didjX)/2. 
Let us focus on the term containing q2 alone. Then wc have 



aq2{t)n[{didjTZ){didjX)]/2 



[-aq2 (d^x) idjn)id^d,n) - aq2 (d^x) n {d^d^n)]/2 

aq2 {(diX) (diTZ) (d^TZ) - TZ didjUdiTl) {djX)]}/2 + aq2Qnil{d^n)/2 . (D5) 



For the term (5^7^) [diX) {S^TV) the following equalities hold 
c(i) {diR) {diX) {d'^n) = -cd^kTld^XdkTZ - cd.TZd.kXdkTZ = 2cnd^kTldikX + cRd.d'^ndiX + cRdi'Jldid'^ X 



so that 



= 2cndikndikX - cdiUdiXd^n - cJld^nd^X + cUdiUdid^X , 

c{t) {diR) {diX) {d'^n) = cUdikTldikX - cQn{d^n)n/2 + cQn{din){din)/2 

= cn{didjn){didjX) + cQn{din){diR) + cQn{dnf /2 , 

where we have defined dijA = didjA. Putting all these partial results together, it follows that 



(D6) 



(D7) 



■ n{^^^Jn){^^^JX) + — (At - 2A5L1) n{didjn){didjX) ^ 



= —a 



2AqLi ^ d f At - 2A^Li 



Wi 



dt 



2w- 



+ H 



A7 - 2AsLi 
2w^ 



2 



Wl 



n{didjn){didjX) 



aq2 



- (At - 2A^L^)n{d,d,Ti){d,d,x) + ^ {{d,x) {d,n) (d^n) - nd,d,[{din) {d,x)]} 



2wi 
a 



A9L1 Qii {dnf . 



2wi 

We choose q2 such that 

2^6^! d / At - 2A5L1 



Wl 

In this case we have 
2A6L1 



dt 



2wi 



H 



At ^ 2A5L1 
2wi 



92 AgLi 

2 Wl 



3H 
2^ 



(At - 2A5L1) . 



(D8) 



(D9) 



^d^UdiTZdiX 



^ n{didjTZ){didjX) + — (At - 2A5L1) TZ{didjn){didjX) ^ 

(At - 2A5ii) [n{didjn){didjX) + 3Hn{didjTZ){didjX)] + ^ {{diX) {d{R,) (d-'n) - ndidjUdiTZ) {djX)]} 

2A^Li) (iZdiU - TZdiiZ) diTZ + TZTZid'^n) - — A9L1 QTZi diTZ) i diTZ) 
2w\ 2 Wl 



-— A9L1 QTZidTZf. 
2wi 

Consider the following combination 



(DIG) 



ac{t)n{d,djn){d^djX) + 3aH c{t)^Z{^^^J^Z){^^^JX) 

= -ac{d.iTZ) [{djTZ){d.,djX)+TZd^d^X] - 3aHc{d^TZ) [{djTZ){didjX) +TZd;d^X] 

= -acXd^dj[{d,TZ){djn)] + ac{d^X) [{dTZf +TZd^TZ] - 3aH cX d,dj [{d,TZ) (djTZ)] +3aHc{d^X) [{dTZf +TZd^TZ] 

= ac{dnf [d^{X + 3HX)]-ac{X + 3HX)didj[{diTZ){djn)]+acnd'^n{d^X)+3aHQcnn^^ 

. c d 



a? dt 



{a^Q 7t) {{dTZ)^ - d'^didj [{d^'JZ){djTZ)]} + acTZd'^TZ {d^X) + 3aH QcTZiZd'^TZ , 



where we have used the properties 

X = Qd-'^']Z, and 



d 



d^{X + 3HX)= (— + 3H] {Qn) = a-^ — {a'^QTZ). 



dt 



d 



dt 



(Dll) 



(D12) 
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The reason why we have introduced this combination is that it is a part of the hnear equation of motion (j29[) . Then 
Eq. ((DTO)) reduces to 



'^^^^'R{d;djn){didjX) + (At - 2A5Ll)il{^^^J'R){^^^jX) - — - 



-( 

•wi 
aq2Q 



aq2 
2 



— (A7-2A5Ll)7^a^7^(9^A') 
2^1 2^1 



Q (Aj - 2A5L1) 7^7^a^7^ - ^r-^ (A7 - 2^5X1) (7^ai7^ - 7^aJ7^) 9,7^ 

2wi 



nn{d^n) AqLi Qn{din){din) AqLi qtz {dny 

2 wi 2wi 



(D13) 



Combining these terms with the other ones in the a/2 term, we find 
a/2 = 



2?i;ia2 dt 



{a'QK) ({OK}' - a-'a,aj[(aR)(3jK)]} + ({dtx) (a,K) (3^k) - KaajKaK) (3^^)]} 



Wi 2?i;i 2 



7^7^a27^ - — Agii Qii{&Rf + a Asn(dn? 

2wi 



+ a(A2- A3L1 + As^^ - Aj^] n'^d^n - (A7 - 2^5X1) n (d^TZ) (d^x) 

\ Wi wi J 2wi 

- ^{Ai- 2A5L1) ('JZ^^n - n^^'R.) diil - — A9L1 Qn{^^'R){^^n). 
2wi wi 

This is similar to the formula (131) of Ref. [63|- Let us define 

A = nizd^n , B = ii{dTzf , c = n [dnf , 
T = ndrUdrn , n = n^^il^^n , x = nd^ndiii , 

Then the following relations hold 

B = ~A-I, V = -2T, /C = -2I, J = -2C, cB=-cC-2cI, 

cE ^ cn{d^n){Qn + Qn) = cQA-n[{cQ + cQ)nd^n + cQnd^n + cQnd^n] 

= -cQA-cQV + cQd.iil {TZdiTZ + TZdiTZ) = Q{-cA -cV + cH + cF) . 
Because of these identities, we can rewrite Eq. (|D14p in the form 



(D14) 



(D15) 



(DI6) 



a/2 = {(^^)' - d-^d,dj[id,n){d,n)]} + ^ mx) {d,n) (d^n) - 7^9,5,[(9i7^) (d^x)]} 

wi 2wi 2 dt\ Iwi J Iwi 

- ^^^n{dTzf + aAsn{dnf + a (A2 - A3L1) n'^d^n - aA9^^7^(a,7^)(a,7^) . (D17) 

2 Wl Wl 

Using the relations Ag = —2wi, cl = —cC + cA, cB = cC — 2c A and the definition of q2, we obtain the form of a/2 
given in Eq. 



By defining 



we have 



Appendix E: Integration of the term a /i 



w = c{t)n{dkiX){dMX), 

Wk = c{t) n {djX){dujX) ~ c{t)dkn{d,x){djX) + c{t) xdj[{djn){dkX)] , 



(El) 
(E2) 



2dkWk -W- 2cndkXdkd^X - 2cXdki{dk'RdiX) + cii{d^kX){d^kX) = ~cll{d^kX){d^kX) . 



(E3) 
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It follows that 

cit) n{d,kX){d,kX) = -cii{d^kX){d,kX) + 2cn{dkX){dkd^x) + 2cXdkiidkndiX) 

= ~cn{d^kX){d^kX) + 2cXdki{dkndiX) - 2c{dkn)idkX)d^x - 2cHid^x){d^x) . (E4) 

There is also another relation 



d{t)nidkiX)idkiX) = -d{dkX)idin)idkiX)-d{dkX)nidkd^x) 

= -dX dki[idkX){din)] + d{dkX){dkn){d^X) + dn{d'^X){d^X) . 



(E5) 



We study the contribution of the last term in Eq. (|60p . i.e. {ciTZ + C2Tl){didjX)(didjX), where ci = A^/wf and 
C2 = Ae/wl. Using Eqs. ((E4| and ([ES]) . we have 

a''cinid,kX)id,kX) + a^ {c2-pi)nidkiX){dkiX)+a^Pin{dkiX)idkiX) 

= -a^{ti+3Hci - pl)^Z{^^kX){^^kX) + a'^[2ciXdki{dkndiX) - 2cidkndkXd^X - 2cin{d^X){d^X)] 

+ a^[{pi - C2)x dki[{dkX){din)] + (c2 - pi){dkX){dkn){d^x) + (c2 - pi)n (d^x){d^x)] , 



where we added and subtracted the term a^piTZ{dkiX){dkiX). Let us define pi such that 

Pi - C2 = 6Hci , 



(E6) 
(E7) 



in which case 



a^'ci'lZid.kXMkX) + (c2 - pi) nidkiX){dkiX) + pinidkiX)idkiX) 
= -a3(ci + SHci - pi)TZ{d^kX){d^kX) 

+ 2cia^{{X + 3HX)dki{dkndiX) - {dkn){dkX)d^{X + 3HX) - n{d^X)[d'^{X + iHX)]} 

= -a^{ci + 3Hci^pi)n{d,kX){d,kX) - 2cij^{a^Qil) {{dk'R){dkX) - d-''d,dj[{d{R){djX)]} 



2ciQnn-^(a^Qn) . 

dt 



where we used Eq. (|D12p . The last term of Eq. ()E8P is expressed as 

dt a-^ dt dt \ a'^ I 



= a^ci Q^il^ + a" (ci - 3Hci)Q'^ 7^7^ 



2 'O'ji 2 



(E8) 



(E9) 



The integrand a^/i can be written as 
a^fi - (Ai+As-^^-qiQ^n^ +a'' 



Ai-Ae^ + Q^- ( ^ 



a 

Wl 



'Ae 


I 




3HA^' 


wf 


dt 




+ 2^ 

wf 



d_ fM 

dt \wl 
2Ar d 



SHA^Q^ 



nid,kX)id,kX) - ^ ^(a^Qn)[{dkTi)idkX) - d-^d,d,[id,n){d,x)]] 



Q 

) — 



ndiUdiX , 



(ElO) 



where we added and subtracted the quantity a^qiQlZ^ . Finally we integrate this term by parts as follows 

d 



a^qiQTT = -n 



qia^Qn^ + gi — (a^Q)7^2 + 2qia^Qnn 
dt 



-nn 



qia^QU + qi — {a^Q)n + 2qia^Qn 
^ d 



-nn 



2qij^{a^Qn) - qij^{a!^Q)n + qia^QTZ 



-2qinnj^{a^Qn) + a^[qi{Q + 3HQ) - Qqi]nn? . 



(Ell) 
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Then Eq. (IeTO)) reads 



3ff.45Q2 



qiiQ + 3HQ)-Qq 



Q 



-a^ ( + A3— - ) - 2q{R.n^{a^Qn) + " 



d f A 



3HA. 



dt 



Wi 

2A5 d 
wf dt 



nd^nd^x 

[a^Qii)[(dkn){dkX) - d-''d^dj[(dji)(djX)]].(Ei2) 



Wi 



The vahie of qi is chosen to match another term in a/2, see Eq. (j7ip . Using the fohowing relation 

c{t) {^jX){^^^Z){^^jX) ^ -c{d^'R){dXf 12 , 
into the following equality 

c{t)TZ{d,kX){d,kX) ^ -c{djX)[d^nd^jX + ndjd^X] ^ c{d^n){dXf /2 + c{d^X)[n{d^X) + {^^n){^^X) 
= c{d^n){dX)^ /2 + cQ^nn^ + cQ'R{d^TZ)id^X) , 

we finally obtain the expression of a^fi given in Eq. ([TO)) . 



(E13) 



(E14) 
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